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• Abstract 
+-> ■ 

In the context of perovskite quantum paraelectrics, we study the effects of 
a quadrupolar interaction J g , in addition to the standard dipolar one J^. 



X 



(3) 

We concentrate here on the nonlinear dielectric response Xp > as the main 



response function sensitive to quadrupolar (in our case antiquadrupolar) in- 
teractions. We employ a 3D quantum four-state lattice model and mean- field 
theory. The results show that inclusion of quadrupolar coupling of moderate 
strength ( J q ~ \Jd) is clearly accompanied by a double change of sign of Xp* 

from negative to positive, near the quantum temperature Tq where the quan- 

(3) 

turn paraelectric behaviour sets in. We fit our Xp to recent experimental 
data for SrTiC>3, where the sign change is identified close to Tq ~ 2>7K. 
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I. INTRODUCTION 



Classical perovskite ferroelectrics such as BaTi03 are very widely studied. As is well 
brought out by recent quantitative ab-initio studies [|l]-|3|], formation of a local dipole moment 
corresponds to formation of a shorter, partially covalent Ti-0 bond while dipole-dipole 
interaction is partly of Coulomb origin (long-range) and partly resulting from the overlap of 
atomic orbitals (short-range). 

The situation is considerably more open and interesting in those cases, mainly SrTiOs 
(also KTaOs) where quantum fluctuations are important to the extent that they actually 
suppress ferroelectricity altogether ("quantum paraelectrics" - Miiller and Burkhard f4|). 
Experimentally, the onset of quantum paraelectric regime in SrTiOs at a "quantum tem- 
perature" Tq ~ 30 -T- 40 K has revealed a variety of intriguing features. They are connected 
in particular with a surprisingly abrupt off-center local displacement of the transition metal 
ion, without ferroelectric or structural long-range order, leading to phase transition-like 
features in EPR f§, EXAFS f§, X-ray and sound velocity A similarly abrupt 



phenomenon has been observed also in NMR studies of KTaOs [11]. In SrTiOs, various 



spectroscopic anomalies have also been found in Raman and Brillouin studies [12 



The reason why quantum fluctuations are so important here is not because the ions 
are particularly light. Rather, it is because the lattice packing is so tight, as compared, 
e.g. to BaTiOs, to leave very little room for the Ti ion to move off-center and form a 
preferential Ti-0 dipole. As the lattice is compressed, the classical ferroelectric Ti off-center 
equilibrium displacement gets smaller and smaller, and the system approaches the classical 
displacive limit, characterized by a vanishing classical Curie temperature T c . In SrTiOs, the 
extrapolated classical off-center Ti displacement and Curie temperature are ~ 0.03A and ~ 
37 K respectively, against ~ O.lA and ~ 400 K in BaTiOs. 111 such a situation, quantum 
fluctuations, even if due to a heavy ion such as Ti, can well remove the ferroelectric long-range 
order and cause a persistence of the paraelectric state down to T — > 0. We have previously 
discussed this scenario in some detail [22[ and also provided a quantum Monte Carlo study 
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of a lattice model illustrating this phenomenon |^TJ . A very recent ab-initio quantum Monte 
Carlo study of SrTiOa and BaTiC>3 has moreover indicated that the exquisite sensitivity of 
the ferroelectric order parameter to this kind of quantum fluctuations is likely to be related 
to the very anisotropic, quasi-one- dimensional electric dipole correlations, which are absent 
for other "neutral" structural order parameters [13 . 



The theoretical study of structural phase transitions and ferroelectricity has so far mostly 
been conducted on models with an anharmonic on-site potential and a bilinear intercell 
interaction (standard model of structural phase transition, fl^JlqQ , as well as on more 



specific shell models [fl6| , pTT|] . Likewise, the studies of quantum paraelectricity have been 



performed within this model [18,19|, or some approximate version of it, like the Ising model 



in transverse field p0| or the quantum 4-state clock model ||21|| . All these studies predict 
the existence of a critical value for the strength of quantum fluctuations, above which static 
ferroelectricity disappears in the ground state. The crossover from the high temperature 
classical paraelectric to low temperature quantum paraelectric is predicted to be a smooth 
crossover, without any sharp features. These models are therefore unable to account for the 
experimentally observed abrupt phenomena at Tq ~ 30 4- 40 K in SrTiCV 

While a thorough understanding of what is going on is still missing, it seems clear that 
on cooling below Tq the individual Ti ions move off-center inside their cell with surprising 
abruptness, whereas their mutual ordering between different cells only progresses slowly as 



T is lowered and never becomes total, down to the lowest temperatures ||22|| . An abrupt 
onset of the local ferroelectric lattice distortion is not well explained by the standard model 
even if the local coupling of the polarization to the lattice strain is taken into account, as 
done for example by Zhong and Vanderbilt In the latter quantum study, no particular 
sharp features have been observed on cooling the system down to T ~ 5 K, contrary to the 
experimental observations. It therefore appears highly probable that in order to achieve a 
satisfactory description of low temperature regime of SrTiC>3, the extension of the standard 
model to include some physically new ingredient is required. 

A very natural extension is, apart from introducing a coupling of the soft mode to other 



phonon modes, to drop the limitation of a simple, dipolar interaction cos(0j — <f>j) and take 
into account also higher angular terms (0j is the phase angle of a 2-component XY order 
parameter in cell i, see section 2). The next order term with square symmetry has an angular 
dependence cos2(0j — (f>j), typical of the quadrupolar interactions. [] 

This paper is devoted to a study of the effects of such a quadrupolar-type of interaction 
on the onset of quantum paraelectricity in SrTiC>3. Deferring to section 2 a more detailed 
presentation of a simplified lattice hamiltonian which takes both dipolar and quadrupolar 
coupling into account, we note that the impact of the new coupling might be of qualitative 
importance. In fact, as we pointed out in quantum effects in a quadrupolar system are 



very different from those in a dipolar system. A purely quadrupolar quantum system exhibits 
reentrance in its phase diagram, while a dipolar one does not. Reentrance in quantum 
quadrupolar systems is beautifully exemplified by the data of Silvera et al. on ordering of 
HD under pressure ||23|| . It also appears very straightforwardly in simple models, at least in 
the mean- field theory p2||2"4"|| . 

The ferroelectric perovskites, unlike HD, are unlikely to be dominated by the quadrupolar 
interaction. Our point, however, will be that they exhibit features which cannot be explained 
by purely dipolar couplings. Experimentally, direct detection of quadrupolar couplings in 
a ferroelectric system is not completely straightforward, since the main quantity, the linear 



dielectric susceptibility Xi; is unaffected by them, at least on the mean-field level p5 
However, the third order nonlinear dielectric susceptibility Xp\ defined by P = Xp^E + 
Xp^E 3 + . . ., turns out to be very sensitive to quadrupolar couplings. In particular Morin 
and Schmitt |25] (in a completely parallel magnetic context) showed that Xp^ wm change 
sign from negative to positive, as a system reaches a temperature region where quadrupolar 
effects become important, while it will stay negative so long as they are unimportant. 



1 This is not the full form of the true Coulomb quadrupole-quadrupole interaction, much as cos((/>j — 
t>j) is not that of a true dipole-dipole interaction. These forms are adopted here for their simplicity. 
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Very recently, new detailed data have been obtained by Hemberger, Lunkenheimer, 
Viana, Bohmer and Loidl pB| , who measured Xp i n SrTiOV They find strong structures 

(3) 

in Xp between 30 and 60 K, precisely in the region near Tq, where the Ti ions suddenly 
move off-center on cooling, giving rise to quantum fluctuating local dipoles, and quantum 
paraelectricity sets in. 

(3) 

In this paper we study this problem, and show that the behaviour of Xp can be under- 
stood if we assume the existence of quadrupolar interactions in SrTi03, besides the standard 
dipolar ones. In section 2, we first introduce a simple quantum lattice model including dipo- 
lar as well as and quadrupolar couplings. In section 3 we calculate Xp within mean-field 
theory, and find that depending on the relative strength of the quadrupolar and dipolar 

(3) 

interactions, Xp exhibits a rich behaviour as a function of temperature. In section 4 we 
consider the application of this calculation to the case of SrTiCV By comparing the cal- 
culated temperature dependence of the third order dielectric susceptibility to the observed 
data |26| we estimate the strength of effective quadrupolar interactions in SrTiCV The 

(3) 

results indicate that in SrTiC>3 a double change of sign of Xp should take place, above and 
below T ~ Tq, with decreasing temperature. The upper one from negative to positive Xp ls 

(3) 

related to onset of quadrupolar effects; the lower one from positive back to negative xp ■> to 
the onset of quantum paraelectricity. Before concluding, in section 5, we shall briefly discuss 
the possible origin of the quadrupolar interactions in connection with the lattice dynamics 
of SrTi0 3 . 



II. SIMPLIFIED LATTICE MODEL OF SRTIO3 

The hamiltonian is quite generally a function of 3 x 5 = 15 continuous coordinates 
per cell in the cubic phase (5 atoms per unit cell) or of even more coordinates in the 
tetragonal antiferrodistorted state. Even simplifying to include only polarization, strain 
and antiferrodistortive degrees of freedom, we are still forced to include at least 3 + 3 + 3 = 9 
degrees of freedom per cell. A proper quantum treatment, for example via PIMC ||T3[ of 
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so many continuous degrees of freedom is at present very demanding. A quantum lattice 
model, although much less realistic, and incapable of describing the displacive behaviour of 
a perovskite in its classical regime (T ^> Tq) is more directly amenable to study either via 
PIMC, or, much more simply, via mean-field theory. 

In this section we introduce a simplified discrete lattice model of a QPE system including 
dipolar as well as quadrupolar interactions. A discussion of these simple models in the 



context of quantum paraelectricity of SrTi03 has already been presented in Ref. |21| , [22 



where a variety of lattice models has been introduced and considered. For convenience we 
briefly repeat here the main ideas. 

We start from the experimentally observed XY-character of the incipient ferroelectricity 
in tetragonal SrTiCV The physical reason why at zero stress and field the polarization is 
confined to the (x, y) plane is given by the presence (below 105 K) of the antiferrodistortive 
order parameter, which expands the lattice in (x, y) plane, while contracting it along the 
z-axis, taken along (001). The onset of ferroelectric polarization is favored in the expanded 
(x, y) plane, and disfavoured in the contracted z-direction. In the (x, y) plane, the polariza- 
tion has four easy directions, namely (±100), (0±10). The corresponding quantum four-state 
clock model introduced first in |2T[ and considered in various versions in [[EJ appears thus 
to be a minimal model for the description of SrTiC>3 at sufficiently low temperatures. The 
model neglects the radial degrees of freedom associated with the continuous displacement of 
the Ti ion from the center of the oxygen cage in the cell % of the crystal. It assumes that the 
displacement is of fixed magnitude and can be completely characterized by a discrete plane 
rotor angular variable (pi. This variable is allowed four possible values (pi = 0, 7r/2, tt, — 7r/2, 
and the corresponding four quantum states labeled as |0), |1), |2), |3) constitute a basis of 
the on-site Hilbert space. The quantum effects are mimicked by allowing the clock variable 
(pi to hop onto its two nearest orientations, i.e. from fa into fa ± | with amplitude —t. This 
is expressed by a hamiltonian term Hi° p , represented in the site basis by the matrix 
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0-10 



1 




-10-10 



(1) 



0-10 



1 



1 



1 



For simplicity, interaction between different cells % ^ j is limited to nearest neighbours, and 
the adjacent cells are supposed to interact via the usual "dipole" interaction 



While the origin of H d is widely discussed, that of H q is not, and this is to our knowledge 
the first place where such a term is invoked for a displacive ferroelectric. The strength of 
the H q term in SrTiC>3 will be estimated in section 4 and its possible origin will be discussed 
in section 5. We note that the classical four-state model H = H d + H q is the well-known 
Ashkin- Teller model [p^j . The 2D classical A-T model is already considerably richer than 
the XY model (^). For the sake of simplicity, we shall treat our model on a 3D cubic lattice, 
and use a mean-field theory. This approximation, while far from accurate, is at least not 
totally unacceptable in (3+l)D and is attractive due to its great simplicity. 

The complete hamiltonian that will be considered in the following is then given by 



H d = -J^cosO;-^), 



(2) 



as well as via a "quadrupole" interaction 



H q 



Jq^COs2(0i - 4>j) . 

(ij) 



(3) 



H = J2H-° P + H d + H q . 



(4) 



III. MEAN FIELD THEORY OF THIRD ORDER DIELECTRIC 



SUSCEPTIBILITY 



In this section we start from our microscopic hamiltonian (^) and derive a free energy 
expansion, using a mean-field approximation. From the free energy expansion we then 
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(3) 

calculate the third order dielectric susceptibility as a function of temperature. Our 



treatment will parallel very closely that of Morin & Schmitt [25] for the magnetic case. 

(3) 

We will show that the behaviour of Xp & t l° w an d intermediate temperatures relative to 
the quantum temperature Tq is profoundly modified by the presence of the quadrupolar 
interactions. 

We take as our mean-field hamiltonian a sum of on-site mean-field hamiltonians, 



# = E^ = E(^-^.-^) 



(5) 



corresponding to a single quantum rotor in external fields e and q, coupling to the on- 
site polarization and quadrupole moment P X i and R4, respectively. For simplicity, we shall 
assume from now on that the polarization component P y is always zero and omit the subscript 
x on P x ; we also set \t\ = 1 as the energy and temperature scale. The matrix of the on-site 
hamiltonian (||) in the basis of the 4 states |1), |2), |3), |4) reads 



(6) 



q -1 -1 
-le-g-1 
-1 q -1 
-1 -1 -e-q 

The trial free energy F t (e, q) per site as a function of the mean- field variational parameters 
e, q is given by 



F t {e,q) 



where 



(H - H°) = F - QJ d l -P 2 - QJ q l -R 2 + eP + qR 



F = ~Tr e"^° , 



(7) 



(8) 



and P, R are now a ferroelectric and a quadrupolar order parameter. The hamiltonian (||) 
cannot be diagonalized analytically, but the expansions of the energy levels in powers of the 
fields e, q around their zero values (we assume that our system, like SrTiC>3 in the absence of 
fields, is paraelectric) can be found by perturbation theory. For our purpose it is sufficient to 
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find the expansion to fourth order in e and to second order in q. The following expressions 
are found for the eigenvalues of 

B » = - 2 -i e2 + ^ 4 -(i + | e V + je 2 ? + ... (9) 
Ei = -q (10) 
E 2 = (l-±e 2 )q + ... (11) 

^ = 2 + i e2 -S e4 + 'i + | e V + J^ + ... (12) 

whence the power expansion of free energy (|3|) can be calculated, 

P — a/W^ 2 -v /r 2 -x/ (2 ^ 2 /r -^ 3 L 4 _L. 



where 



Xo 1} = ^tanh/3 

1 + 4/3e 2 ^ - e 4 ^ 3 
Xo =■ 



X2 



16(1 + e 2 ' 3 ) 2 
1 - - 
2(l + e 2 P) 2 



x &> = -Itaah* . (14) 

The polarization and quadrupolar order parameters P and R can now be calculated as 
derivatives of fllTf ) with respect to the fields e and g 

P = -^ = X^ + 2 X ^g + X ^ 3 + ... (15) 
R=-^ = X2q + X ( 2 ] e 2 + ... . (16) 

In order to obtain an expansion of the free energy term Fq in terms of the order parameters 
P, R, we must invert the expansions ([15]) and (|Tj]) and express the fields e, q in terms of 
powers of P and R. We find the expressions 



X?' (rf>) 2 *2 (xf') 4 V * 

1 ( 2 ) 

g=— i?-— P 2 + ... , (18) 

X2 (x£ ] ) 2 X2 



which after substituting to (p~3|) yields the desired form 



r'l 



(19) 



At this point, we can express the full free energy F t (|7|) as an expansion in powers of the 
order parameters, namely polarization P and quadrupole moment R. It reads 

1 1 v (2) 

Ft(P,R) = ^(xSV - QMP 2 + ^(xi)- 1 - QJ q ]R 2 - -ih—P 2 R 

1 1 (Xo ) X2 

As in real SrTiOa at zero stress and field, we assume P = R = 0. In the neighbourhood of 
this minimum of F t , we now wish to eliminate R, yielding an effective free energy expansion 
in powers of the polarization P alone. Minimizing ( p0|) with respect to R, we obtain 

(2) 

R = m X2 P 2 > ( 21 ) 

and after substituting this expression back to fl20"| ) we obtain the expansion of F t in terms 
of P only 

Ft(P) = i [(xf »,-. _ WJP . _ ( x f + H^gj p. . (22) 

Now it is worth noticing that the coefficient of P A has been renormalized downwards (i.e. to- 
wards negative values) by the quadrupolar coupling J q . This follows because Xo 3 '' is negative 
definite, while the other term dependent on J q is positive as long as J q < (6x 2 ) _1 - 

By adding an interaction term —EP, corresponding to an external electric field E, to 
the free energy fl22|) and minimizing with respect to the polarization P, we can calculate 
the first order (linear) and third order (nonlinear) dielectric susceptibilities Xp an d Xp\ 
defined as usual by the expansion 

P = X { p ) E + X { p ) E :i + ... . (23) 

In our model we find 
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1 - 6 J dX { o ] 



xf = - " or- U 3) + ™tH ■ (25) 



The nonlinear susceptibility (|25|) can also be written in a slightly more transparent form 

X?=X?(J, = 0) + ^f^, (26) 

where we made use of the relation (x^) 2 = (Xo 1 "*) 4 ) verified by Eq. (|T4|) . While the first 
order dielectric susceptibility Xp\ m the mean field approximation, is independent of J q , 

(3) 

we see that the third order susceptibility Xp is a sum °f its original value in the absence 
of quadrupolar coupling which is always negative, and of an additional positive term. The 

(3) 

latter can eventually overbalance the negative term and reverse the sign of xp ■> if Jq is 
large enough. The temperature dependences of Xp f° r Jd = 0.273 and a series of values of 
quadrupolar coupling J q are plotted on Fig.l. We see that a) the effect of renormalization 
of Xp^ is most strongly pronounced at intermediate temperatures T ~ Tq, and b) there is 

(3) 

a region of values of J q where Xp now turns positive at intermediate temperatures, while 
staying negative at sufficiently high and low temperatures. Finally, if J q is stronger, there is 

(3) 

a temperature where xp * 00 an d an independent (anti) quadrupolar ordering takes place. 

(3) 

We point out here that positive Xp > which means a negative overall coefficient in front 
of P 4 term in (P2|) would imply the ferroelectric transition at the classical level to be of 
first order. The QPE state for positive Xp can therefore be seen as a frustrated first order 
transition. 



IV. APPLICATION TO THE CASE OF SRTIO3 

We now wish to apply the above results to the case of SrTiC>3. As mentioned in the 
introduction, SrTiC>3 belongs rather to the class of displacive than to the order-disorder 
ferroelectrics. The use of a continuous model should thus be preferred to a discrete one, if 
a realistic comparison to experiment is to be attempted. However, the main purpose of this 



(3) 

paper is to pursue the qualitative differences induced in the behaviour of Xp by a presence 
of a quadrupolar interaction with its different angular dependence, and these features should 
be correctly reproduced also by a discrete model, which in turn has the big advantage of 
being easily analytically tractable on a simple mean-field level. It is clearly beyond the 
possibilities of such a discrete model to describe the sharp dipole onset observed at Tq ~37 
K, since the local off-center displacement is described just by the radial part of the local 
mode which is left out completely in the present discrete model. 



Before attempting any comparison between experimental data [26 and our results (|24l) 
and (p5|), we have to convert the latter formulas from dimensionless units back to normal 
ones. We consider first the linear dielectric susceptibility. If the displacement of the Ti ion 
from the center of the cage is associated with a dipole moment /z, the dipole density per 
unit volume being n, then x~p can be conventionally written in form of the familiar Barrett 
formula 



X { P = T T B T m > ( 27 ) 

P Q coth £L - T n ' y J 



derived in many different contexts [^9|,[|]. The temperatures To and T\ are given by To = 
3Jd/kB,Ti = 2t/ks, and have the meaning of classical Curie temperature and of quantum 
temperature, respectively. This formula has been used to fit the Xp* data on SrTiOs 0,|26 



and produces a reasonable fit over a temperature range from to 100 K with T\ ~ 88 K, 
T ~ 36 K, and ^fV ~ 1q5 K - The 

microscopic hamiltonian parameters corresponding to 
these values of To and T\ are Jd = 12 K, t = 44 K. We note in particular that the ratio 
Jd/t = 0.273 is not far from the critical mean-field value of (Jd/t) c = 1/3, where our model 
undergoes a transition from para (J < J c ) to ferro (J > J c ). 

2 

A problem arises, however, with the prefactor ^ , since with one dipole per cell n = 
1/a 3 = 1.68 x 10 28 m~ 3 , where a = 3.9 A is the lattice constant of SrTiOs, the required value 
of the dipole moment per cell is \i — 3.4 x 10~ 29 Cm. If we take the effective charge to 



be Z* ~ 8 Pm, the corresponding displacement should be d ~ 0.26A, which is an order of 



magnitude larger than the value of ~ 0.03A realistically estimated in Ref. fU, as well as 
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in ab-initio calculations ||. The discrepancy is, in our opinion, related to the quantitative 
inadequacy of the discrete model to SrTiC>3, as well as to the mean- field treatment. The 
Barrett formula reproduces the essential features of the temperature behaviour of the linear 
susceptibility well enough, but fails at the quantitative level. 

Let us now come to the nonlinear susceptibility, and the related quadrupolar effects. The 

(3) 

complete expression we have just obtained for \p is 

w/J 4 f / 107/ (2)\2\ 

x f = Lm + -^TT 1 ) • (28) 

(t-6J dX i>>)' \ t-6J, X , ) 

/ON 

where (3 = t/(kpT). We can now try to fit this formula to the Xp data measured on 
SrTiOa ||26|| . For this purpose, however, the form in which the data are presented in Ref. 
26| , namely loglxp^l vs. temperature plot is not so well suited, since the mere knowledge 



of the magnitude of the complex quantity \Xp \ without its phase angle, mixes together the 

real and imaginary part and conceals the desired information about the sign of the real 

(3) 

part. Our static mean-field theory produces, of course, a real static Xp ■ We have therefore 



attempted to extract Re x^p from the data [26|, making use of the corresponding phase 
angles kindly provided by J. Hemberger [31]. The resulting experimentally derived Re x$ 



(corresponding to the lowest electric field intensity E = 50 V/mm) is shown on Fig. 2 as a 
function of temperature. While at low temperatures this quantity is negative, on heating it 
crosses zero at T = 33 K where it becomes positive, passes through a maximum at T = 45 
K and then approaches zero again, from the positive values. Above T ~ 60 K, it probably 
turns small and negative again, although experimental uncertainty appears larger in this 



regime. With hindsight, we re-examined also some old data of Fleury and Worlock |32| for 

(3) 

SrTiC>3 and found that evidence for a positive Xp at T = 40 K can be evinced there, too 
(Fig. 5 b, Ref. H). 

Having in mind the difficulties related to the prefactor, encountered in case of Xp 1 we 

(3) 

decided to allow this to be a free parameter. Anticipating that the best fit of Xp does 
not necessarily have to lead to the same values of Jd and t as the fit of Xp\ we decided to 
preserve the value of Jd/t = 0.273 and allow for J q and t to be free parameters. The result of 
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this fitting is the curve on Fig. 2, which corresponds to t = 75 K, Jd = 0.273t = 20.5 K and 

(3) 

J q = 4.8 K. With these values, the positive maximum of Xp is well fitted, although at lower 
temperatures the fit is less perfect within the present theory. At higher temperatures on the 
other hand it is hard to assess the agreement between the theory and data; in any case the 

(3) 

theory predicts \p to become negative again at high temperatures and approach zero from 
below. From Fig. 2 it is clear that by assuming a quadrupolar interaction of strength about 
5 K, the theory can qualitatively account for the temperature dependent double change of 
sign of Xp ■ Quantitatively, the agreement is worse than in case of Xp\ which is perhaps 
not surprising for a thermodynamic quantity which is a higher order derivative of the free 
energy. The required value of t is now larger by ~ 70 % and the value of the displacement 
required in order to obtain the necessary value of the prefactor is even an order of magnitude 
larger than in case of Xp\ namely 2.3 A, which is clearly unphysical. As in the case of Xp > 
we attribute these discrepancies to the crudeness of the model. On the other hand, a firm 
result of this calculation and fit is that the presence of J q appears to be crucial in order to 
account for a positive x^'p °J intermediate temperatures. 



V. POSSIBLE ORIGIN OF THE QUADRUPOLAR INTERACTION IN 

PEROVSKITES 

Before closing, we discuss the possible microscopic origin of the effective (anti)- 
quadrupolar interaction J q in perovskites. Recently, we did in fact invoke for SrTiC>3 |22[ 
the possible importance of an effective quadrupolar coupling J q between the dipoles in dif- 
ferent cells, in addition to the dipolar one Jd- The sign of this interaction J q is negative 
("antiquadrupolar"), such that dipoles prefer to be either parallel or antiparallel, to being 



orthogonal. One physical source anticipated in Ref. could be related to a lattice strain 
effect, as follows. When a dipole bond appears in one cell, it will cause an elongation of that 
cell along the dipole direction, and a shrinking in the orthogonal direction via electrostric- 
tion. This deformation, to avoid strain gradients, will extend to the neighbouring cells where 
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it will favor equally well a parallel or an antiparallel dipole orientation, but will disfavor the 
orthogonal one. The coupling then arises so as to reduce strain gradients associated with 
orthogonal dipoles, but absent for parallel or antiparellel dipoles. In order to estimate the 
strength of the resulting interaction, these intuitive considerations have to be formalized. 

Formally, the term representing the coupling between the local polarization and strain 
is, for a perovskite structure, linear in strains and quadratic in polarizations, and according 
to Ref. |§ it reads 

H mt = \ E E B la p Vl (Ri)M^)M^) , (29) 

1 i la/3 

while the elastic energy is given by 

H el = WY, C ikVi(Ri)Vk(Ri),. (30) 

i Ik 

— * — * 

where r)i(Ri) and u a {Ri) are the strain and local soft mode amplitude components at site i, 
respectively, Bi a p are the strain-ferroelectric soft mode coupling constants and are the 
elastic constants. The total strain at site i consists of the homogeneous component arising 
from the uniform deformation of the whole system as well as of the inhomogeneous local 
strain. Both components have to be taken into account in a different way, because the 6 
local strain components per cell are not in fact independent quantities, but rather suitable 
linear combinations of the 3 independent acoustic displacement components per cell. To 
proceed, one has to express the local strains in both ( ^9|) and (0) in terms of the acoustic 
displacements. Integrating out the acoustic displacements from the total hamiltonian, one 
finds the effective interaction between the polarizations in different cells, which turns out to 
be fourth order in polarizations and contains also the quadrupolar interaction which we are 
interested in. In our case, however, such procedure would not be as straightforward as in the 
corresponding classical case, because the complete hamiltonian contains the quantum kinetic 
energy term which does not commute with the coupling term containing the polarization. 
The simpler classical calculation should, however, provide at least a rough estimate of the 
strength of the effective J q induced in this way. 
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Before doing this calculation in detail, however, it is actually possible to estimate an up- 
per limit to the strength of the effective J q by means of a simple argument. The induced in- 
ter site interaction term of the form (u(Ri) ■ u(Rj)) 2 = \u{Ri) 2 u{Rj) 2 {l + cos2{(j)i — (j)j)),i ^ j 
certainly must be weaker than an analogous induced term for i = j, (u(Ri)-u(Ri)) 2 = u A {Ri), 
which represents a renormalization of the on-site quartic term u A (Ri). This renormalization 

B 2 

is clearly of the order of In order to get an upper estimate for the interaction strength 
J q (fl), we must also multiply by |<i 4 , where d is the typical value of the local displacement. 
Taking the typical values of B\ a p ~ 1.4 hartree/bohr 2 and Cik ~ 5 hartree (from Ref. ||), 
together with the displacement d ~ 0.03A, we arrive at a value about 0.5 K, which is an 
order of magnitude lower than the value of 5 K estimated in the previous section to be 
necessary to explain Xp i n SrTiC>3. 

A more accurate classical calculation of J^ tram is most conveniently performed in the 
Fourier representation, due to the translational invariance of the hamiltonian. The long- 
wavelength limit of the problem has been studied in the 70 's when it was shown that the 
resulting interaction is nonanalytical in fc-space and has a long-range tail in the real space 
5B|jr4]] . A similar 2D calculation of indirect strain-induced coupling has been performed for 



the case of SrTiC>3 [33], taking into account all the Fourier components from the Brillouin 
zone, which allows extraction of the induced interaction for each pair of lattice sites i ^ j. 
We will not go here into the full detail of that calculation and merely quote the main result, 
which is that the strength of the induced interaction is of the order of 0.1 K, a factor of five 
lower than the previous upper limit estimate. We conclude that the origin of the effective 
quadrupolar interactions of the required strength is not likely to be elastic, and should be 
sought elsewhere. 

Another possible source could be long-range interactions. Coupling one cell with further 
neighbours will generally yield terms with higher angular Fourier components. The soft 
phonon mode dispersion away from k = is in fact very non-sinusoidal, precisely due to 
Coulomb interactions |35|| . A third possibility might finally arise from anharmonic terms 



like (u(Ri) ■ u(Rj)) 2 of non-strain origin. At this stage we are not in a position to decide 
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among these possibilities. It would therefore be highly desirable to extend ab-initio studies 
similar to those in Ref. || beyond the lowest order considered there ("local anharmonicity 
approximation"), and try to determine also the required higher order terms of the energy 
expansion. 

Because of the close relationship between these expansions and the harmonic and anhar- 
monic phonon properties, which in case of KTaC>3 and KNbC>3 are also quite well reproduced 
by extended shell models |36|], it might also be possible to use them for this purpose. 



VI. CONCLUSIONS 

(3) 

In conclusion, we have considered the nonlinear dielectric susceptibility xp ■> an d its 
behaviour with temperature for a model system exhibiting quantum paraelectricity below a 
"quantum temperature" Tq. The model is chosen so as to include a quadrupolar coupling in 
addition to the normal dipolar one. It is found that the effect of the quadrupolar couplings 

(3) 

is to cause a temporary switch of sign of Xp from negative, well above and well below Tq, to 
positive at T ~ Tq. Recent data for x^p °f SrTi0 3 appear to agree well with this predicted 
behaviour implying an effective quadrupolar coupling J q ~ \ Jd- Work which remains to be 
done includes 

a) better understanding of the true microscopic origin of J q , presently still unclear 

b) understanding the possible relationship of this new coupling to the exciting and largely 
unexplained phenomenology observed near Tq, in Miiller's original experiments ||, and those 

(3) 

which followed it. For the time being, we have noted that a positive Xp at Tq implies that 
the QPE state results due to quantum frustration of an otherwise first-order ferroelectric 
transition. 
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FIGURES 

(3) 

Fig.l. Third order nonlinear dielectric susceptibility %p as a function of temperature 
for a fixed value of dipolar coupling Jd and various values of quadrupolar coupling J q . 



Fig. 2. Fitting of the experimentally measured Re(xp^) [0)0)0] °^ SrTiO"3 with formula 
f). Full points correspond to experimental data, solid line to the theoretical curve. 
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Re(x (3) P ) [(V/mm)" 2 ] 




